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Abstract 

Using a high statistics sample of photo-produced charm particles from the FOCUS 
experiment at Fermilab, we report on the measurement of the ratio of semileptonic 
rates r(D+^r M +^) = 625 ± .045 ± 0.034. Allowing for the Kn S-wave interfer- 

ence measured in [1], we extract the vector to pseudoscalar ratio J^ ^^ _o - " '-— — 



=o 



r(D+-iTV + &v) 



0.594±0.043±0.033 and the ratio S^^ = 1.019±0.076±0.065. Our results 
show a lower ratio for yW^kI T than has been reported recently and indicate the 
current world average branching fractions for the decays D + —> K (/t/ + , e + )u fjL)e are 
low. Using the world average B(D + — > K^t: + t: + ) [2] we extract B(Z) + — > K fi + is) = 
(9.27 ±0.69 ±0.59 ±0.61) %. 
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**Corresponding author; present address: Vanderbilt University, Dept. of Physics 
and Astronomy, 6301 Stevenson Ctr., Nashville TN, 37235 (USA). Tel.: 615-343- 
8295 (Cell: 303-579-6032). E-mail: will.johns@vanderbilt.edu 



1 Introduction 



There is currently some controversy concerning the relative rates of the charm 
vector semileptonic decays that proceed via a K* and the charm pseudoscalar 
semileptonic decays that proceed via a K. Theorists originally expected the 
rates to be about the same [3], but other theoretical predictions and experi- 
mental measurements in the 90's tend to favor a smaller vector semileptonic 
rate (see Tables 1 and 2). 

A more recent experimental result [4] indicates that the ratio of rates is closer 
to unity than measured previously. Rather than measure the rate for the 
pseudoscalar and vector semileptonic decays directly, as was done in previous 
measurements by the same experiment [5,6], the result uses the average from 
the PDG for the determination of B(D + — ► K e + u e ) [2]. Other experiments 
measure the semileptonic rates for D + and D° decays and form a ratio of 
vector to pseudoscalar rates (see Table 1). 

Since the D + and the D° particles only differ by the light quark, exclusive 
semileptonic rates for the decays of these particles are expected to be equal 
through SU(3) symmetry. A comparison using the current world averages of 
the pseudoscalar decay branching fractions along with the D + and D° lifetimes 
[2] indicates that the pseudoscalar semi-electronic rates (the error on the D + 
pseudoscalar semi-muonic rate is too large for a meaningful comparison) are 
different at the 99% confidence level: T(D + -► K°e + u e ) -T(D° -»■ K~e+v e ) = 
—25 ± 9.7 ns _1 . This result is surprising and merits further investigation. We 
intend to show in this paper that the difference in rates and the recent CLEO 
[4] result are in part due to the pseudoscalar semileptonic branching fraction 
reported in the PDG [2] for the decay D + — > K°£ + U£ being too low. 

We measure the ratio — — ~^_ o^ — directly, using decays with a very similar 

topology. Previous measurements of this ratio used comparisons between the 
D + and the D°, relied on PDG branching ratios, and/or used decays with 
different topologies (for instance where one of the modes requires an added 
pion). By reconstructing the neutral kaon in the microvertex detector of the 
FOCUS experiment [7,8] through the decay Kg — > 7r + 7r~, we take advantage 
of the studies and work performed to produce precise lifetime measurements of 
the long-lived charm particles. By measuring the D + decay, we take advantage 
of the extensive work performed to understand the decay D + — > K~7T + fi + u tl 
[1,9,10] and our result is the first to incorporate the interference effects de- 
scribed in [9] and measured in [1]. 

The data for this analysis were collected using the Wideband photoproduction 
experiment FOCUS during the 1996-1997 fixed-target run at Fermilab. The 
FOCUS detector is a large aperture, fixed-target spectrometer with excellent 



vertexing and particle identification used to measure the interactions of high 
energy photons on a segmented BeO target. The FOCUS beamline [11] and 
detector [7,8,9,12] have been described elsewhere. 

Table 1 

Previous results compared to the FOCUS result. Notice that some results [13,14,15] 
are admixtures of different charm species related through the isospin argument, 
[14,15] are correlated since the same D + — > K [i + v result is used, and in the 
CLEO(02) result [4], the T{D+ -» K°e + iy)/T Total comes from the PDGOO [16] esti- 
mate for T(D+ -» ifVz/)/r Tota i. 



Experiment 


Quantity 


Result 


CLEO(91)[5] 


r(D°-*K*-e+i>) 

F(D u ^K-e+u) 


0.51 ±0.18 ±0.06 


CLEO(93)[6] 


T(DQ->K-e+v) 


0.60 ±0.09 ±0.07 


CLEO(93)[6] 


V{D+~,K* y 'e+u) 


0.65 ±0.09 ±0.10 


E691(89)[13] 


F(D u ^K-e+u) 


0.55 ±0.14 


E687(93)[14] 


r(D+->x*V + ^) 

F(D°^K- l i+u) 


0.59 ±0.10 ±0.13 


E687(95)[15] 


T(D+-+K*"n+v) 

T(D ^K-[i+v) 


0.62 ±0.07 ±0.09 


CLEO(02)[4] 




0.99 ± 0.06 ± 0.07 ± 0.06 (±0.12) a 


FOCUS(04) 


T(D+-+K*"n+v) 


0.594 ± 0.043 ± 0.030 



a The PDGOO [16] error for T(D+ -► K°£ + iy)/T Total , omitted [18] in the CLEO [4] 
result, is shown in parentheses. 



2 Event Selection 



We identify D +, s through the 3-body decay D + — ► h~n + n+v^ (where the h 
represents a pion or a kaon and charge conjugate modes are implied throughout 
this paper). To search for candidate events, a search strategy common to both 
modes is employed. This is possible since roughly 10% of all Kg's reconstructed 
in the ir + 7i~ channel are reconstructed using hits from the FOCUS silicon 
microvertex detectors. 1 Where possible, we have chosen cuts similar to those 
optimized in other FOCUS analyses to enhance charm content and particle 
identification. 

Two opposite sign tracks reconstructed using information from the silicon 
detectors are required to form a vertex with a confidence level exceeding 1%. To 
suppress background from short-lived, primarily non-charm particles produced 



1 The technique used in [17] to reconstruct D + — > Kgir + using the bulk of the Kg 7 s 
is not applicable in this case due to the missing neutrino. 



Table 2 

Theoretical predictions for the Vector to Pseudoscalar ratio. No distinction is made 
for electrons, muons or charm species. In some cases, listed with a (*), the ratio was 
calculated from information given in the publication. 



Model 


T{D^K*l+v) 
F(D-^K£+u) 


WSB(85*) [3] 


1.13 


LFR(90*) [19] 


0.7,0.67 


SUMR(91)[20] 


0.50 ±0.15 


LAT(92*) [21] 


0.86 ±0.27 


LAGR(93*) [22] 


0.56 


LAT(94) [24] 


1.1 ±0.6 ±0.3 


MITBAG(94) [23] 


0.56 


ISGW2(95) [25] 


0.54 


QPT(96) [26] 


0.65 


RQM(96) [27] 


0.57 


QM(97) [28] 


0.62 


LFM(97) [29] 


0.68 


SR(97*) [30] 


n 47+ - 19 


QM(00),DISP(01) [31,34] 


0.63 


COVQ(01*) [32] 


1.01,0.72 


CVLFD(Ol) [33] 


0.66,0.64,0.67 


PCL(Ol) [35] 


0.54,0.63,0.57,0.67 


EFT(02) [36] 


0.5 ±0.2 



in the targets, this 2-track vertex is required to occur at least 1 standard 
deviation outside of target material. The two tracks are then formed into a 
single track that is combined with a candidate muon to make a putative D + 
vertex with a confidence level exceeding 1%. This D + vertex is required to 
occur at least 1 standard deviation outside of target material as well. 

Due to the relatively long lifetime and Lorentz boost of charm candidates, 
the primary interaction vertex and secondary D + decay vertex can have a 
significant separation along the beam direction. Tracks not used to form a D + 
are used to construct a set of candidate primary vertices with confidence level 
greater than 1%. The highest multiplicity primary vertex candidate with the 
highest separation, £, from the secondary in units of error, ag, greater than 3 is 
retained. Since the significance of the separation between the interaction and 



decay vertices is an essential tool used to separate charm from background, 
we perform a scan in this variable between £jat = 3 and £/<j£ = 23 (roughly 
a factor of 3 in yield) to judge the stability of our result. We find that our 
results become stable after a cut of £jci > 5 is reached. Additionally, the fit 
quality for the determination of the D + — > Kg^v^ yield becomes optimal 
(x 2 /DOF < 1) between £/<j£ = 11 and £/<?£ = 19. For the final sample, a 
value of £/<j£ = 13 is used. The ratios of the yields of individual modes to 
that expected from simulations exhibit the same stability. This is a strong 
indication that our cuts are effective at removing any short-lived (non-charm) 
backgrounds that could mimic signal. 

Muon candidates are required to be within the acceptance of the inner muon 
system in FOCUS [9]. Since the rate of muon misidentification increases at 
low momentum, we require that the momentum of muon candidate tracks be 
greater than 13 GeV/c. Muon candidates are required to have associated hits 
in the muon system sufficient to meet a minimum confidence level of 5% for 
the muon hypothesis where at least 5 of 6 planes of the detector must record 
hits consistent with the candidate track. 

To separate decays proceeding through the K^tt+jj+v^ channel, from those 
decaying through the KgH+v^ channel, we use additional vertex, particle ID 
and invariant mass requirements. 

The vertex representing the K~tt + is required to occur within 3 standard 
deviations of the (K~7i + )fi + vertex, and a three track vertex formed from the 
K~ ,7r + and fi + tracks must exceed a confidence level of 5%. Since the K s 
lifetime is long compared even to the charm lifetime, we find a very effective 
cut to reduce non-K s contamination is to impose a requirement that the K s 
vertex have a large separation from the charm decay vertex. We require that 
the vertex representing the K s — ► 7r~7r + decay occur at least 15 standard 
deviations downstream of the (7r~"7r + )/i + vertex. 

We use the Cerenkov system [12] to identify pions and kaons. For each track, 
^obs = —2 log (£) is computed, where C is the likelihood that a track is con- 
sistent with a given particle hypothesis. For the track in the K~7i + fi + vertex 
with charge opposite to the muon, we require W ^ s {ti) — W Q ^ S {K) (kaonicity) 
be greater than 2.0, and for the track same charge as the muon in the K~n + fi + 
vertex, we require W \, S (K) — W / b s (vr) (pionicity) be greater than 0.0. For pions 
that form a candidate Kg, we require that the pion likelihood be favored over 
the particle hypothesis that forms the minimum likelihood W min — W i, s (ir) be 
greater than -5.0 (a loose cut). 

Background from D + — »■ K~7t + tt + and D + — ► K s tt + , where a pion is misiden- 
tified as a muon, is reduced by requiring that the visible mass M(h~7r + n + ) < 
1.8 GeV/c . In order to suppress background from D* + — > D°n + — > (K~ ji + u^)n' 



we require M(K~ir + fi + ) — M(K~~ /i + ) > 0.2. In order to enrich the Kg sam- 
ple we require the K s invariant mass be within 2 standard deviations of the 
nominal reconstructed value. 

To reduce backgrounds from topologies consistent with extra tracks com- 
ing from the secondary vertex (such as expected backgrounds like D° — ► 
K 7r - ;U + z/ M or high multiplicity charm and non-charm backgrounds not present 
in the simulation) we reject candidates where any track not used to create the 
secondary vertex can be included in the secondary vertex with a vertex fit 
confidence level exceeding 10%. 



3 Analysis 



We fit the K s fi + invariant mass and the K~n + invariant mass. Fit components 
are a combination of Monte Carlo histograms: one generated with the mode 
of interest and others representing background components. A maximum like- 
lihood is used where the predicted number of events in a bin^ is described 
by: 



7V(predicted)j = PsSignalj + P^Background 



■i j 



where Ps and the P/s (more than one background shape is used) are fit pa- 
rameters. The number of predicted signal events in the data is then described 
by J2f ms PsSignalj where Signal^ is the number of reconstucted Monte Carlo 
events for the mode of interest in a given bin. 

We find that the fit to determine the yield of the D + — > K~7i + fi + h , fJi events 
requires only 2 components: one for signal and one for background. The D + — ► 
K^7i + fi + u ll matrix element has only recently been fully measured, and we 
simulate the signal shape and efficiency with the results from [1]. We represent 
the background shape in the K~rr + mass by generating a Monte Carlo in 
which we simulate all known charm decay backgrounds while removing D + — ► 
K~'K + fj, + u fl from the generated particle mix. 

We have tried several approaches to fitting the invariant Kgfi + mass his- 
togram. To generate D + — > K /i + f^ events, we use the simple pole form for 
the pseudoscalar to pseudoscalar semileptonic form factor where: 



f±(q 2 ) = /±(0)/(l - q 2 /M 2 pole ), q 2 = (P D - P kaon ) 2 = M^_ virtual . 
We use Mp i e = 2.11 GeV/c and /-//+ = —0.7. Since we accept almost the 



entire Kgfi + invariant mass range, small changes in the choice of M po i e and 
/_//+ have a negligible effect on our final result. To represent the background 
in the data, several techniques were investigated. 

If the background in the data is primarily due to D — *■ K n^v^ when a 7r° 
or 7r + is missed and D + — > K /i + f M when either the Kg or the fi + is not 
from the D + , it should be sufficient to perform a fit including only these two 
components. We find that this approach produces consistent results but poor 
fit quality at low £ fat's. Fit quality improves though at the cost of stability 
if the 2 (or more) lowest mass bins are removed from the fit. We also find a 
slight improvement in fit quality and stability at high £/<j£ if a Breit-Wigner 
compor 
the fit. 



component centered at 892 MeV/c with a width of 50 MeV/c is added to 



In order to try and improve the fit quality, we added a background shape 
to the previous 3 component fit by generating a Monte Carlo in which we 
simulate all known charm decay backgrounds while removing D + — > Kg^v^ 
from the generated particle mix. We found that the fit quality did not improve, 
the error on the returned fit increased, and the results were not stable below 
£/a = 9. Even though the results agreed quite well at higher £/a, we decided 
to investigate this behavior further. 

In order to reproduce a background shape which may contain unsimulated 
backgrounds, we took events which had at least one extra track consistent 
with the secondary vertex at vertex fit confidence levels between 30% and 
90%. In order to gauge the specific effect of non-simulated backgrounds on 
the final result, invariant mass histograms were formed from both the data 
and from a Monte Carlo in which we simulate all known charm decays except 
D + — ► K /U + z/ M and D + — > K 7r°/i + z/ M . Using these shapes in addition to 
the 3 component fit significantly improves the fit quality. The x 2 /DOF for 
the fit using the background from the data is acceptable (~ 1 or less) at 
all £/a£S. There is now a component of the signal in the data background 
subtraction. As a consequence, the signal and background are correlated in 
the fit, and the fit errors increase. In the fit using the simulated background, 
we find that the x 2 /DOF is about 1 unit higher than the data represented 
background fit until i/cr^s above 11, where the fit quality becomes equivalent 
between the two representations. The difference in x 2 /DOF at low £/a/s is 
likely due to non-charm or short-lived charm decays, appearing primarily at 
low Kgfi mass according to the results of the binning tests, which are not 
included in the simulation. Even though we see a difference between the data 
and simulated background in x 2 /DOF at low i/crg, the results of the fits are 
in good agreement and stable above Ijot = 3. 

Our quoted result uses the fit to the K s fi + invariant mass with the simulated 
background from higher multiplicity secondary vertices at £/<j£ = 13. It is 



important to note that in all four of the fits described, the 3 component 
fit, the fit with the inclusive simulated background, and the fits with the 
background from higher multiplicity secondaries produce equivalent values 
above i/crg = 11. This is likely due to the background being dominated by 
D — ► K-Kfi^u^ decays at higher £/&£. 

The fit to the data for both modes is presented in Fig. 1. We find 555 ± 
39 Kgfi + u^ decays and 9871 ± 127 K~it + ji+v^ decays. 
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Fig. 1. The fit to the Kg/j, + invariant mass (left) and the K~ir + invariant mass 
(right). The fit to the data (error bars) is shown as a solid line, and the background 
described in the text in shown as a dotted line. For both plots, the difference be- 
tween the solid and dotted lines represents the signal. Note the good fit to the K^fi 
invariant mass and the cleanliness of the K~ir + invariant mass. 



Thus, our ratio — — ^_^ M becomes: 



r^ 



Wn+ Vll ) 



1/2 #^7rV + ^( FIT ) <K° sf i + u, 



2/3 #ATg/x+^(FIT) e(K-K+pi+v, 



-FT*0 



The 2/3 accounts for the probability that K u decays to K"ir + (see above) 
and the 1/2 accounts for the probability that K° decays to Kg. The Kg — > 
7t + tt~ branching fraction is accounted for in the Monte Carlo generation. The 
number of events determined from the fit to the data for each mode is labeled 
(FIT), and the reconstruction efficiency for each mode determined using the 
Monte Carlo indicated by an e. In order to quote the K* component as a 
separate result, we separate the resonant and non-resonant components using 
the technique outlined in [1]. Fully 95.0 ±0.5% of the K~n + sample is proceeds 
through a K . Thus we find: 



T(D+ -► Knu+v) 



T(D+ -»• TTVi/) 
and, 



T(D+ -► i^*V + ^) 



0.625 ± 0.045 



_ 0.594 ±0.043. 

Our systematic tests of the result are outlined in the next section. 

4 Systematic Checks 



Our systematic uncertainty comes from known quantities which are not in- 
cluded in the fit (such as the S-wave contribution estimate), unanticipated 
variations in the data not accounted for in the simulation, and variations due 
to the fitting technique. 

To determine the amount of D + — ► K /tx + z/ M contained in our D + — ► K~ir + ji + u 
signal, we corrected our estimated yield of D + — ► K~7T + fi + v by 0.950 ± 0.005. 
We compute this fraction by integrating over the K~7r + fi + v phase space the 
model intensity and parameters from reference [1] with the S-wave amplitude 
set to zero and dividing this value by the same with the S-wave amplitude and 
phase set to the measured values. The uncertainty (±0.005) is determined by 
varying the amplitude and phase by the errors indicated in [1] and noting the 
difference. 

In order to assess a systematic uncertainty due to larger variations of M po i e , 
we varied M po i e between 1.86 and 2.31 GeV/c in the Monte Carlo generation. 
The resultant simulated signal histograms are then used to repeat the fit 
used to obtain the final result. The sample variance from the returned fit 
results is retained as the systematic uncertainty due to M polc . To estimate the 
systematic uncertainty due to different values of /_//+, we repeated the fit at 
Mp i e = 2.11 GeV/c with /_//+ = 0.7 and kept the difference as the error 
estimate. 

In order to assess a systematic error to the measured ratio from unanticipated 
variations in the data not accounted for in the simulation, we placed a vari- 
ety of pertinent cuts on the data and computed a sample variance from the 
returned values. The cuts described below are in addition to those previously 
applied. Unless a particular meson or mode is mentioned, cuts are applied to 
both modes used to calculate the ratio simultaneously. 
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To check for non-i^ 7r + 7r~ backgrounds, we set the normalized Kg mass cut 
to values of 3 and 1. To check for backgrounds where a neutral long-lived 
particle such as a A is misidentified as a Kg, we made a cut requiring that the 
difference in the magnitude of the K s candidate pion momenta be no greater 
than 70% of the sum. We also increased the Cerenkov requirement on both 
pions so that W bs(K) — W o bs0r) (pionicity) be greater than 0.0. Since this 
latter cut primarily removes low momentum K s 's, it is an effective tool at 
examining the Kg acceptance as well. 

To check for unanticipated backgrounds and differences between the simu- 
lation and the signal mode due to event topology, we investigated how the 
ratio behaves with a variety of requirements on the detailed location of the 
reconstructed event in the spectrometer. We required the D + vertex be lo- 
cated downstream of the first interaction target, downstream of the second 
interaction target, upstream of a trigger counter located near the main silicon 
tracking system, or upstream of a location roughly between the 2 downstream 
target silicon system planes. We also required the Kg vertex be 2 cm down- 
stream of the D + vertex, upstream of the main (last 12 planes) silicon system, 
or downstream of the target (1st 4 planes) silicon system [8]. 

To look for short lived or non-charm background, we increased the requirement 
to 3 standard deviations that the D + vertex occur outside of target material, 
and we required P v isibie > 30 GeV/c . As a check for higher multiplicity decays 
feeding into the signal, we specified that the maximum allowable confidence 
level that an additional track be consistent with the secondary be 1%. 

In order to specifically look for decays feeding into the signal where a particle 
is misidentified as a muon, we chose cuts that should reduce the probability 
of contamination while leaving high efficiency for signal. We placed a cut 
on the muon requiring that the momentum measured by both magnets in 
the spectrometer be consistent. We increased the muon momentum cut to 
> 20 GeV/c . We also increased the requirement on the muon identification 
confidence level to 15%. 

To check for additional background in the Kn mass, we increased the Cerenkov 
likelihood difference cut on the kaon from > 2.0 to > 4.0, and increased the 
cut on the mass difference used to cut out D*'s to 0.25 and 0.30. 

We assess a systematic uncertainty from these cut tests by computing a sample 
variance of the returned values for — — ^_„ M — . The larger contributions to 
this systematic error estimate are listed in Table 3. 

In order to test the fit, we performed repeated tests where we Poisson fluctu- 
ated the bins of both the data and fit histograms and performed repeated fits. 
Our tests indicate that the K s ji + u^ yield error is underestimated by 6% (due 
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to finite Monte Carlo statistics), and the statistical error in the final result is 
boosted to account for this difference. We also looked at the standard devia- 
tion between the four different fits tried. Although the fits agree remarkably 
well at I/up = 13, we felt a more conservative approach was to choose a value 
that was common to four returned ratios below and above the chosen £/ag. 
The systematic uncertainty estimated from the four fits is added in quadrature 
to the previously described estimates to assess a total systematic uncertainty 
in the ratio (see Table 3). 

Table 3 

Larger contributions to the systematic uncertainty in the ratio — - — ^_a - — . Note 

that the uncertainty associated with the non-S-wave Kir contribution applies only 
to the estimate of the K* fraction in the Kir signal. 

Systematic Contribution Value 

Normalized K Q S Mass Cut 0.008 

Secondary Vertex Location 0.017 

K° s Vertex Location 0.013 

Muon Magnet Consistency 0.012 

Muon Momentum Cut 0.008 

Total contributions from cut variations 0.028 

Mp i c and /_//+ variation 0.015 

Contribution from fit variations 0.013 

Total systematic uncertainty 0.034 

S-wave Fraction (K* ratio only) 0.003 



5 Summary and Conclusions 



Our result represents a substantial improvement over previous results for the 
ratio of the vector to pseudoscalar decay in the muon channel. We find 



T(D + -» Knfi + u) 



T(D+ ^ICfi+is 



0.625 ± 0.045 ± 0.034 



and 



T(D + -► K*°fi+u) 



T(D+ -> TCfi+v 



0.594 ±0.043 ±0.033. 
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Where the first error is statistical and the second error is the systematic uncer- 
tainty. The K* J K result agrees with older results, but disagrees with a recent 
result from the CLEO collaboration (with no corrections (~ 3% effect) made 
for phase space and form factors). One possible source of this difference can 
be due to the vector decay as is discussed in Reference [1]. Another source is 
detailed below. 

By using our vector /pseudoscalar result and our measurement of the branch- 
ing ratio r ( D +~^ K -J+J+\ [10], correcting for the updated values of the S-wave 
contribution [1], we determine: 



T(D+ -> ICu+is) 

1 - 1.019 ±0.076 ±0.065, 



T(D+ -> K~n+n+) 

where we have added the statistical and systematic errors separately in quadra- 
ture. 

Using the PDG [2] values for the absolute branching fraction of the decay 
D+ -► K-71+71+: (9.1 ± 06)%, we calculate, 



B(L> + -► K% + u) = (9.27 ± 0.69 ± 0.59 ± 0.62) %. 

The third error is due to uncertainty in the D + — > K~7t + tt + branching fraction. 
Our result is a substantial improvement over the present world average [2] of 
7.0i|o%. Besides the difference between the D + and D° semi-electronic rates 
mentioned in the introduction, there are other reasons to believe that the PDG 
value for the D + semi-electronic mode 6.5 ± 0.9% is low. 

We can compare the D° and D + semi-muonic rates as we did for the semi- 
electronic rates in the introduction. We find good agreement with T(D + — ► 



TC ii^Vfj) — T(D° — ► K /tx + z/ M ) = 11 ± 11 ns x (where no correction for the 
difference in phase space between D + and D° has been made). 

It is likely that the semi-muonic rate is lower than the semi-electronic rate by a 
few percent, and this is consistent with what is measured for the isospin conju- 
gate decay of the D° using values from [2]. One sees r L ~^._ h^ryrW^K^^) = 
1.12 ±0.07, for the D° but ^^^^^^ = 0.72 ± 0.11 for the 

D + using the result in this paper. The PDG estimates that this ratio should be 
around 1.03[2]. This is reasonable since a very large positive ratio for /_//+, 
which increases the semi-muonic rate, is not likely given the E687 [15] re- 
sult for /-//+, and radiative corrections for the semi-electronic mode, which 
can lower the measured semi-electronic rate, are expected to be small. It is 
also unlikely that any such large, unanticipated corrections apply to the D + 



13 



semi-electronic mode exclusively. 

Finally, it is interesting to note that the sum of the rates measured by FOCUS, 
corrected as suggested in the PDG [2] to estimate the semi-electronic modes, 
is B(D+ -► ((l.05)Kn + (1.03)^°)^+^) = 14.9 ± 1.2%. This is closer to the 
current world average inclusive electronic rate D + — »■ e + (anything)= 17.2 ± 
1.9% than B(D + -> ((3/2)^-^+ + K )^^) = 12.9±{;|% (both from [2]), 
suggesting that a large, previously unseen, semileptonic decay mode for the 
D + is unlikely. 
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